RKKY interactions in the regime of strong localization 
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We study the influence of strong nonmagnetic disorder on the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactions between diluted magnetic moments in metals. We find that the probability 
distribution for the RKKY interactions assumes strongly non-Gaussian form featuring long tails. 
Since such distributions cannot be characterized by its moments, we define a typical value of the 
interaction amplitude, which we find to be exponentially suppressed in presence of Anderson local- 
ization. Our results present a plausible and physically transparent picture describing how Anderson 
localization effectively eliminates the long range nature of the RKKY interactions. 
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I. INTRODUCTION 

It has long been appreciated that locaHzed magnetic 
moments in metals interact through indirect RKKY in- 
teractions mediated by the conduction electrons .A In a 
clean metal, the RKKY interaction has a long-range 
oscillatory part, with an amplitude which decreases as 
a power law of the distance between the impurities, 
I{R) ^ cos{2kFR)/R'^ {d is the dimensionality of the 
system and kp the Fermi wave vector). This behavior 
is well understood to be a direct consequence of the ex- 
istence of a sharp Fermi surface characterizing itinerant 
electrons. 

In presence of impurities and disorder, this behavior 
may be substantially modified. It is well known that suf- 
ficiently strong disorder can lead to multiple-scattering 
processes which can trap the electrons through the pro- 
cesses of Anderson localization. In this regime, one may 
expect that the long-range character of RKKY interac- 
tions should be suppressed, refiecting the reduced mobil- 
ity of conduction electrons. The essential physical ques- 
tion is how this process precisely takes place as the dis- 
order strength is gradually increased and the system ap- 
proaches the localized regime. 

The infiuence of weak nonmagnetic disorder on the 
RKKY interactions has been studied in considerable 
detail,— >^>^>^>^ and is by now very well understood. These 
studies have established that the main effect of weak 
disorder is to randomly modify the phase of the RKKY 
oscillations due to impurity-induced phase shifts of the 
electronic wave functions. As a result, the RKKY in- 
teraction decreases exponentially when averaged over 
disorder, {I{R)) e~^/', where I is the mean free 
path. Early work^ thus predicted that the range of 
RKKY interactions becomes essentially cut off at R^ I, 
and can be neglected at larger distances. More careful 
considerationi^iS' discovered that this naive argument is 
incorrect, essentially because most relevant quantities do 
not depend on the average value of the interaction, but 
instead on the typical value of its amplitude. Since this 



quantity is not affected by the random phase shifts, it 
is essentially unaffected by weak disorder, and decreases 
with the distance R in the same power-law fashion as in 
the clean system. Therefore, the presence of weak dis- 
order is not expected to lead to any significant changes 
in the physical properties which are dominated by the 
RKKY interactions. 

The infiuence of stronger nonmagnetic disorder on 
RKKY interactions has been, however, so far poorly ex- 
plored. In an important study, Lerner examined the 
probability distribution of the RKKY interactions in the 
metallic phase in the presence of strong disorder within 
generalized nonlinear a model and 2 + e expansion J It 
was found that the quantum interference corrections do 
not change the power-law decay of all the even moments 
of the interaction distribution, which remains the same as 
in the pure metal, but make the coefficients attached to 
these moments increase critically with disorder. As the 
Anderson localization regime is approached, the higher 
moments increase much faster than the variance, which 
therefore no longer represents a good characterization 
of the typical interaction strength. Although this work 
clearly points out to the importance of strong fiuctua- 
tions of the RKKY interactions, it does not explain how 
does the distribution function evolve as one enters the 
regime of Anderson localization. 

Within an Anderson insulator the electrons are bound 
to impurities, and thus can hardly be expected to gener- 
ate the long-range part of the RKKY interaction. How 
can we have at the same time large moments of the dis- 
tribution of RKKY interactions, along with its fast de- 
cay with the distance? At first glance these two argu- 
ments seem inconsistent and the situation confusing and 
paradoxical. Clarifying these issues is an interesting and 
important problem, since one expects the magnetic cor- 
relations to play a crucial role in the physics of disorder- 
driven metal-insulator transitions in general. 

The resolution of this puzzle is in fact quite simple, 
as we explain in this paper. We find that in presence of 
strong disorder and localization, the distribution func- 
tion develops a strongly non-Gaussian form featuring 
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long tails. In such cases it is well known that all mo- 
ments of the distribution can assume very large values, 
while at the same time a typical width of the distribution 
can remain very small. Instead of the arithmetic aver- 
age (i.e. standard deviation of the distribution) we find 
that the typical value of the interaction is better char- 
acterized by the geometric average of the distribution, 



ItypiR) = 



^ p(iln[/(i?)]^> 



This quantity is exponentially 



suppressed in the presence of Anderson localization, ex- 
plaining how the long-range part of the RKKY interac- 
tion is suppressed. 

To illustrate these ideas and obtain quantitative and 
reliable results, we numerically study the distribution 
of the RKKY interactions within the Anderson insula- 
tor phase. Although our numerical results are obtained 
within one dimensional model, we argue that the same 
concept of the typical value of the RKKY interaction 
can be used to physically explain the qualitative change 
of the form of RKKY interactions from long ranged to 
short ranged during the disorder driven metal-insulator 
transition in general dimensions. In the rest of this paper 
we present a detailed description of the model, followed 
by the numerical results and discussion. 



II. RKKY INTERACTIONS IN A DISORDERED 
METAL 

The interaction energy between two local moments Si 
and S2 embedded in the metallic host at ri and r2, re- 
spectively, is given by the Hamiltonian 



H,nt = -J'Si •S2x(ri,r2), 



(2.1) 



where x is the zero frequency nonlocal electronic suscepti- 
bility and J is the exchange coupling constant. Therefore 
the calculation of the interaction between diluted mag- 
netic impurities reduces to the calculation of the suscep- 
tibility which, expressed through the Matsubara Green's 
functions, is given by 



2 X - 

X(ri,r2) = -2^G'a,„(ri,r2)G^„(r2,ri) 



(2.2) 



Here w„ is the fermionic Matsubara frequency, and f3 is 
the inverse temperature. 

In the clean case, i.e. in the absence of the nonmagnetic 
disorder, in three dimensional (3d) system, and for R 
much larger than the lattice spacing, the susceptibility is 
equal to 



XoiR) = - 



2mkp cos(2fci?i?) 
(27r)3i?3 



(2.3) 



where i? = |ri — r2| and m is the effective mass. 

In the presence of weak disorder the phase of xiR) 
becomes random and x{R) averaged over the disorder 
configurations is exponentially suppressed 



(X(i?)) = XoiR)e-^^\ 



(2.4) 



where I is the mean free pathf^i^iii^i^ The second mo- 
ment (variance) of the probability distribution, however, 
remains long ranged and has the same power law depen- 
dence as in the clean system^i^i^ 
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(2.5) 



Jagannathan et al^ have found that t he sq uare root of 
the forth moment of the distribution, y/ (x^), is compara- 
ble in magnitude to the second moment. Therefore, the 
susceptibility distribution is non-Gaussian, but its typi- 
cal value is well characterized by the square root of its 
second moment. The same remains true in 2d as welli^ 

This statement is, however, not valid in the regime of 
strong disorder. As shown by Lernerj^ using the general- 
ized nonlinear a model and performing a 2-l-e expansion, 
further increase of the disorder results in very rapid in- 
crease of all the even cumulants of the distribution. More 
precisely, for the cumulant ((x")) of the order n 



((X"(i?))) 



R 
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(2.6) 



where the parameter u > 1 as the disorder is increased 
and the system approaches the Anderson transition.™ In 
this case the typical value of the distribution cannot be 
determined by the value of its moments. 



III. NUMERICAL RESULTS 

In order to examine the form of the RKKY interac- 
tions in the regime of strong disorder, we proceed to a 
numerical study. We consider a tight binding model 



H = -t 



E 

{«i>o" 



cLc 



'JO- 



h.c.) + gjC, 



(3.1) 



with nearest neighbor hopping t and on-site random po- 
tential Ei, which is distributed uniformly in the interval 
[~W/2,W/2]. We calculate the interaction between the 
magnetic impurities embedded into this system at dis- 
tance R. From Eqs. i|2.ip and (|2.2p we see that all the 
information that we need is contained in the single parti- 
cle Green's functions G'a;„ (ri, r2). In the matrix notation 
the Green function 



H) 



(3.2) 



and the problem reduces to the numerical summation of 
the corresponding matrix elements over Matsubara fre- 
quencies LOn which will be done in the zero temperature 
limit. In order to obtain good statistics with large num- 
ber of disorder realizations and to reduce the finite size 
effects, we concentrate in the following to a one dimen- 
sional system. 

The numerical results for the clean system. Fig. 1, re- 
produce the well known oscillatory form of the electronic 
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Figure 1: The average susceptibility as a function of distance 
R in the presence of weak disorder of strength W. The chem- 
ical potential is /i = —0.96. 
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Figure 2: Scaled average susceptibility R{xiR))- With the 
l/R dependence gone, we can see that disorder introduces a 
damping factor to the average interaction strength. 



susceptibility with the power law decay with the distance, 
x{R) ^ cos{2kpR)/R. Here the power law exponent is 
equal to 1 since we are working within one dimensional 
(Id) model. (The disorder strength W and the chemical 
potential fi are measured in units of half the bandwidth 
2t, and the distance R in units of the lattice spacing. Our 
system had 500 lattice sites.) We then consider the sus- 
ceptibility (x(-R)) in presence of weak disorder averaged 
over hundreds of disorder configurations. The average 
susceptibility weakens as the disorder is increased. 

In Fig. 2 we remove the l/R dependence by multiply- 
ing the average value (x) by R. For weak disorder {x)R 
follows an exponential decay as predicted long time ago 
by de Gennes.— However, what we are really interested 
in are the probability distributions of the electronic sus- 
ceptibility in the presence of stronger disorder. 

Fig. 3 shows the probability distribution of the scaled 



Figure 3: Probability distribution P{xR) for W = 0.35 and 
fi = —0.96. The distribution width is clearly dependent on 
R. 

susceptibility P{xR) in the presence of strong disorder, 
and for several values of R. As in the remaining part 
of the paper, the results are obtained by averaging over 
hundreds of disorder configurations on the lattice with 
500 — 1000 lattice sites. The width of the distributions 
are dependent on R. More importantly, cursory exam- 
ination of the data indicates that the distributions are 
distinctly non-Gaussian. This feature is important for 
the following analysis, since the non-Gaussian shape of 
the distribution prevents us from using the standard de- 
viation as a measure of the the width of the distributions. 
Instead, we define a typical value of the width as the ge- 
ometrical average 

X*,p(i?)^e<^>">^^(«)). (3.3) 

Fig. 4 shows a comparison of standard deviation and Xtyp 
for describing the width of the same distribution. We can 
see that the standard deviation is infiuenced by the long 
tails, and is far too large to be a useful description of the 
width. Xtyp, however, gives a good estimate of the width 
of the distribution. 

Examining xtyp for different strengths of disorder 
(Figs. 5 and 6), we find that for each value of the dis- 
order W 

i?Xt,p(i?) ~e"^/« (3.4) 

for sufficiently large R, where ^ defines the localization 
length. Therefore, if we define an adjusted susceptibility, 

XA(i?) = i?e«/«x(i?) (3.5) 

the distributions of xa{R) for large R will be entirely 
independent of R within each disorder strength. Com- 
paring Figs. 3 and 7 shows how the use of this adjusted 
susceptibility causes these distributions to collapse to a 
single scaling function. We then combine the data for 
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Figure 4: Distribution -P(x) (full black line) for R — 100 and 
W = 0.2. Gaussian (blue dotted line) is taken to have the 
same standard deviation as P{x)- Xtyv (I'^d dashed line) is 
a better measure of the distribution width than its standard 
deviation. 
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Figure 5: 7? Xtyp on linear axes, 
decrease with increasing R. 
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Figure 6: Rxtyv on semilogarithmic axes. Rxtyp features an 
exponential decay with R. 



0.02 



< 



0.01 



1 


\\ — R = 200 




- - R = 250 




1|- ■ ■ R = 300 






^^^^^^^ 






,1 



-20 



-10 



10 



20 



Figure 7: Scaled P{xa) for W = 0.350. After scaling, the 
distributions collapse. 



several distances R, and thereby obtain more precise dis- 
tributions for each value of disorder W. For each disorder 
strength there now is a single characteristic distribution 
(independent on R) as shown in Fig. 8. Interestingly, 
the probability distributions are quite asymmetric in the 
presence of stronger disorder, ferromagnetic interactions 
being much more probable than the antiferromagnetic 
ones. In the strongly localize regime we expect in fact 
the interactions to be ferromagnetic for R smaller or of 
the order of the localization length. The obtained distri- 
butions show this behavior even for R much larger that 
the localization length. 

We then plot the tails of the distributions and find 
that for very strong disorder they become very long, as 
shown in Fig. 9. The tails of the distribution appear to 
converge to a universal power law form in the limit of 
strong disorder. The form of the tails is qualitatively the 



same for positive and negative side of the distribution. 
The existence of such long tails indicates that the width 
of the distributions for higher disorder cannot be accu- 
rately characterized by their moments. The moments of 
such a distribution are extremely large, while the typical 
value is, in fact, very small. Therefore, the long range 
part of the RKKY interactions is strongly suppressed in 
the strongly localization regime. 

We have also compared our result for the mean free 
path I, obtained from Eq. (|2.4p . and the localization 
length ^, from Eq. (|3.4p . Plotting I and ^ as a function of 
W, see Fig. 10, we find that I and ^ are both proportional 
to and find that ^/l « 3.685, which is in a good 

agreement with the analytical result ^/l ~ 44 This anal- 
ysis further confirms the consistency of our interpretation 
of as determined from the decay of the typical RKKY 
interaction amplitude, with the localization length of the 
electronic system. 
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Figure 8: P{xa) for several values of W . The distributions 
are strongly peaked for higher disorder. 




Figure 9: Tails of distributions on log-log axes. For higher 
disorder, the tails become increasingly long. 



IV. CONCLUSION 

In this paper we have examined how the distribution 
function for RKKY interactions becomes modified due 
to Anderson locaUzation effects. We demonstrated that 
the essential effect of locahzation is to exponentially sup- 
press the typical amplitude of RKKY interactions on dis- 
tances longer then the locahzation length, in agreement 
with intuitive expectations. The distribution, neverthe- 
less, remains "broad" in the sense that it develops long 
tails which dominate the statistics. Our numerical re- 
sults thus confirm the analytical predictions of Lerner 
that all even moments of the distribution diverge within 
an Anderson insulator. 

Our results portray an interesting physical picture with 
potentially far-reaching consequences. In the metallic 
regime the RKKY interactions remain long-ranged even 
in presence of weak disorder, and thus a given magnetic 
moment effectively interacts with many others. In this 
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Figure 10: The mean free path I and the localization length ^ 
as a function of the disorder strength. Both I and ^ are found 
to be proportional to . 



regime one may expect a well developed collective be- 
havior of the spin system, leading to magnetic order- 
ing at low temperature. When the Anderson-localized 
regime is approached, the situation is quite different. The 
RKKY interaction between a typical pair of distant spins 
is now significantly suppressed or even neghgibly small. 
Very occasionally, a pair of distant spins will interact 
strongly, due to rare disorder configurations producing 
long tails in the distribution function. If the resulting 
RKKY interaction is antiferromagnetic, then such a pair 
can be expected to lock in a tightly bound singlet - 
thus forming an essentially inert molecule that practi- 
cally detaches from the rest of the spin system. Such 
process are precisely what one expects within the ran- 
dom singlet phase a^i^°i^^i^^ which feature quantum Grif- 
fiths phase anomalies li^. Alternatively, a ferromagnetic 
interaction will lead to the formation of a bound triplet 
state, essentially a spin S = 1 magnetic moment. It is 
interesting to note a degree of asymmetry of the distri- 
bution of the RKKY interactions, that we have found at 
strong disorder. This finding seems to indicate that fer- 
romagnetic correlations may effectively compete with the 
tendency for singlet formation, possibly leading to nano- 
scale ferromagnetismM coexisting with a random singlet 
phase. 

We thus anticipate that Anderson localization pro- 
cesses generically destabilize spin glass ordering in disor- 
dered metals, which is instead replaced by an appropriate 
quantum Griffiths phase. Precisely how these processes 
take place in the vicinity of reahstic metal-insulator tran- 
sitions remains a fascinating open direction for future 
study. 
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